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Abstract 

The distributions that regulate the spatial domains of the Poissonian Voronoi 
Diagrams are discussed adopting the sum of gamma variate of argument two. 
The distributions that arise from the product and quotient of two gamma vari- 
ates of argument two are also derived. Three examples of non Poissonian seeds 
for the Voronoi Diagrams are discussed. The developed algorithm allows the 
simulation of an aggregation of methanol and water. 
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1. Introduction 

A great number of natural phenomena are described by Poissonian Voronoi 
diagrams , we cite some of them : lattices in quantum field theory [H ; con- 
ductivity and percolation in granular composites [1, 0]; modeling growth of 
metal clusters on amorphous substrates the statistical mechanics of simple 
glass forming systems in 2D Q; modeling of material interface evolution in 
grain growth of polycrystalline materials [6[ . We now outline some applications 
in Chemistry. Detailed computer simulation results of several static and dy- 
namical properties of water, obtained by using a realistic potential model were 
analyzed by evaluating the volume distributions of Voronoi polyhedron as well 
as angular and radial distributions of molecular clusters Q. The local struc- 
ture of three hydrogen bonded liquids comprising clusters of markedly different 
topology: water, methanol, and HF are investigated by analyzing the proper- 
ties of the Voronoi polyhedron of the molecules in configurations obtained from 
Monte Carlo computer simulations The local lateral structure of dimyris- 
toylphosphatidylcholine cholesterol mixed membranes of different compositions 
has been investigated on the basis of the Voronoi polygons Molecular dy- 
namics simulation of a linear soft polymer has been performed and the free 
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volume properties of the system have been analyzed in detail in terms of the 
Voronoi polyhedron of the monomers [loj . The molecular dynamics simulation 
of the aqueous solutions of urea of different concentrations are modeled by the 
method of Voronoi polyhedron . 

From the point of view of the theory the only known analytical result for the 
Poissonian Voronoi diagrams , in the following V~P after the two memories 12, 

, is the size distribution in ID . More precisely the theoretical distribution 
of segments in ID is the gamma variate distribution of argument two , in the 
following GV2 IJ]. The area in 2D and the volume in 3-D were conjectured to 
follow the sum of two and three GV2 IJ] . The possibility that the segments , 
area , and volume ofV — P can be modeled by a unique formula parametrized 
with d which represents the considered dimension d{d = 1,2,3) has been recently 
analyzed [l5[. Section [2] of this paper reviews the known formulas of the sum of 
GV2 and express them according to the chosen dimension. Section [3] explores 
the product and quotient of two GV2. Section [4] analyzes three examples of 
Voronoi Diagrams with correlated seeds or non Poissonian processes. Section [5] 
reports the simulation of an aggregation of methanol and water. 



2. Sum of a Gamma variate 

The starting point is the probability density function ( in the following PDF) 
in length , s , of a segment in a random fragmentation 

p{s) — Xe:xjp {—Xs)ds , (1) 

where A is the hazard rate of the exponential distribution. Given the fact that 
the sum , w , of two exponential distributions has a PDF 

p{u) = \^uexp{—Xu)du . (2) 

The PDF of the ID V — P segments , /, ( the midpoint of the sum of two 
segments) can be found in the previous formula inserting u — 21 

p(0 = 2A/exp(-2A0rf(2A/) . (3) 

When transformed in normalized units x — j the following PDF is obtained 

p{x) ^ 2xexj>{-2x)d{2x) . (4) 

When this result is expressed as a gamma variate we obtain the PDF (for- 
mula (5) in Q) 

^^(a;;c) = ^^(cx)'="^exp(-ca;) , (5) 

where < x < oo , c >0 and r(c) is the gamma function with argument c; in 
the case oi ID V ~ P c — 2. As an example Figure [1] reports the histogram of 
length of the normalized Voronoi segments in ID. 
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Figure 1: Histogram (step-diagram) of the V — P normalized segment distribution in ID 
with a superposition of the gamma- variate as represented by equation the number of 
Poissonian seeds is 20000 , d = 1.02 , c = 2.04 , NBIN = 40 and x^=32.52 ,random seeds 



The Kiang PDF has a mean of 



and variance 



1 



(6) 



(7) 



It was conjectured that the area in 2D and the volumes in 3D of the V — P 
may be approximated as the sum of two and three GV2 . Due to the fact that 
the sum of n independent gamma variates with shape parameter Ci is a gamma 
variate with shape parameter c = '^i' area and volumes are supposed to 
follow a gamma variate with argument 4 and 6 3, 13] ■ This hypothesis was 
later named "Kiang conjecture", and equation ([5]) was used as a fitting function 



ISl . Iig . |20|, |2l[, or as a hypothesis to accept or to reject using the standard 

The PDF ([5]) can be gen eralized 

1,2,3)11, Hi 



procedures of the data analysis, see [22, ]2: 



by introducing the dimension of the considered space, d{d 

2d 



H{x;d) 



T{2d) 



{2dx)^'^-^ ejip{-2dx) 



(8) 



Two other PDFs give interesting results in the operation of fit of the V — P 
area/volume. The first is the generalized gamma PDF with three parameters 
(a,6,c) , [2i,[ll,[2l , 



f(x\b,c,d) = c 



exp {—bx'^) 



(9) 
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The generalized gamma PDF has the mean of 



M= ^A^^ ^ , (10) 



and variance 

^ 6-(+r(2±^)r(f)-(r(i±^))^ 



(r(f))' 



(11) 



By the method of maximum Ukehhood the estimators , a, b and c are the solution 
to the simultaneous equations 

-ln(b) - n—^ + V Infx,) = 
c c ^-^ 

i=l 



a 1 c 
n-- - > 



(12) 



n a 
— -n- 

c c 



n—ln{b) + n — fo^J x1ln{xi) — , 



where is the digamma function , n the number of observations in a sample 
and Xi is an observed value. The second one is a PDF of the type (15j 

FiV(x;d) = Consi X cc^exp (-(3d + l)a;/2) , (13) 

where 

^ , V273dTT 

Const = ——J , (14) 

2 23/2 (3 d + r (3/2 d + 1/2) 

and d{d — 1, 2, 3) represents the dimension of the considered space. We will call 
the previously reported function the Ferenc-Neda PDF which has the mean of 

11^1 , (15) 

and variance 

Which distribution produces the best fit of the area of the irregular polygons 
and the volume of the irregular Polyhedron ? In order to answer this question 
we fitted the sample of the area and volume with the three distributions here 
considered . The 2D V ~ P are reported in Figure [5] , Table [T] reports the 
results and Table [5] reports the parameters of generalized gamma deduced in 
two papers with the addition of our values. 

Figure [3] reports the various PDFs here adopted for the normalized area 
distribution in 2D. 
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Figure 2: The V — P in 2D. The selected region comprises 44 random seeds marked by a 
point. 



Table 1: Values of for the cells normalized area-distribution function in 2D. Here we have 
25087 Poissonian seeds and 40 intervals in the histogram; u denotes the degrees of freedom. 



PDF 


V 




k{x]d) {Eq.^) ,d= 1.77 


39 


83.48 


{Eq.m) 


39 


71.83 


G{x-a,b,c) {Eq.^) 


36 


58.9 


h{x) (Eq.l^) 


39 


12636 


r{x) (Eq.^) 


39 


19986 
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Table 2: Generalized gamma parameters a, b and c of the normalized area distribution in 2D 
in different papers with Poissonian seeds. 



reference a b c 

Ferenc- Neda 2007 2.29 3.01 1.08 

Tanemura 1988 and 2003 3.31 3.04 1.078 

this paper 3.15 2.72 1.13 




Figure 3: Poissonian seeds. Plot of the cells normalized area-distribution function in 2D when 
5 PDFs are adopted : k{x ;d) (Eq. ^) ,d=1.87 (full line) , f(x;b,c,d) with data of Ferenc- 
Neda 2007 (dashed), f(x;b,c,d) with data of Tanemura 1988 and 2003 (dot-dash-dot-dash) , 
f(x;b,c,d) with our data (dotted) and f(x;d) (Eq. JTSj) ,d=2, (dash-dot-dot-dot) . 
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Figure 4: Histogram (step-diagram) of the V — P normalized volume distribution in 3D with a 
superposition of the gamma- variate as represented by equation l(8}: the number of Poissonian 
seeds is 1392 , d=3, NBIN=40 and x^=62.63 ,random seeds 

In 313 Figure S] reports the histogram of the volume distribution as well as 
the Kiang's PDF [8] when d = 3 , Table [3] reports the results and Table [4] 
reports the parameters of generalized gamma in two papers with the addition 
of our values. 

Table |4] reports the parameters of generalized gamma PDF deduced in two 
papers with the addition of our values. 

Figure [5] reports the various PDFs here adopted for the normalized volume 
distribution in 3-D. 



Table 3: Values of for distribution of normalized volumes of 3-D cells. Here 21378 Poissonian 
seeds were generated and 40 interval were used in the histogram. 



PDF 




h' 


x' 


k{x;d) (-Eg.®) ,d = 


2.76 


39 


93.86 


f{x;d) {Eq.^,d^ 


3) 


39 


134.15 


G{x;a,b,c) (Eq.l^) 




36 


58.59 
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Table 4: Generalized gamma parameters a , b and c of the normalized volume distribution in 
3D in different papers with Possonian seeds . 



reference a b c 

Ferenc- Neda 2007 3.24 3.24 1.26 

Tanemura 1988 and 2003 4.80 4.06 1.16 

this paper 4.68 3.87 1.18 




Figure 5: Poissonian seeds. Plot of the cells normalized volume-distribution function in 3-D 
when 5 PDFs are adopted : k(x ;d) (Eq. ^) ,d = 2.79 or c = 5.58 (full Une) , f{x;b,c,d) with 
data of Ferenc-Neda 2007 (dashed), f(x;b,c,d) with data of Tanemura 1988 and 2003 (dot-dash- 
dot-dash) , f(x;b,c,d) with our data (dotted) and f(x;d) (Eq.l ll3t ) ,d=3, (dash-dot-dot-dot) 
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3. Product and Quotient of a Gamma variate 

This Section explores the product XY and the quotient X/Y when X and Y 
are GV2. 

3.1. Product 

We recaU that if X is a random variable of the continuous type with PDF , 
f{x), which is defined and positive on the interval < a; < oo and similarly if 
y is a random variable of the continuous type with PDF g{y) which is defined 
and positive < y < oo , the PDF oiV ^ XY is 

Kv)^ r 9{-)f{x)-dx . (17) 

Jq X X 

Here the case of equal limits of integration will be explored , when this is not 
true difficulties arise 2^[2^ . When f{x) and g{y) are GV2 the PDF is 

h{v)= — dx = Z2vKo{A^/^) , (18) 

JQ ^ 

where Ki,{z) is the modified Bessel function of the second kind 27, 2^ with v 
representing the order, in our case 0. The previous integral can be solved with 
the substitution x = | and using the integral representation in [2§| , pag. 183 

1 1 1 

K-^i^) - 2^r^' I 'M---^)^dr . (19) 

The mean of the new PDF, h{v), as represented by formula (fT5]l is 

<v>= w X 32 vKo (4 v^) = 1 . (20) 
Jq 

The mode , m , is at v =0.15067 and Table [T] reports the of the fit of the 
V — P normalized area-distribution in 2D. 

3.2. The quotient 

The PDF oi V = X/Y when X , represented by g(x) , and Y , represented 
by f(y) , are GV2 is 

r{v) = / \w\g(vx)f(x)dx = . (21) 

Jo (l + v) {l + 2v + v^) 

The mean of the new PDF, r{v), as represented by formula (|2ip is 

<v V X 5 dv = 2 . (22) 

Jo {l + vf{l + 2v + v^) 

The mode , m , is at w = and Table [T] reports the of the fit of the Voronoi 
cell's normalized area-distribution in 2D. 
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Table 5: Normalized area distribution in 2D of 166 correlated seeds generated according to 
formula l|23|l . The of data fit when the number of classes is 40 for three PDF 



PDF 


name 






H{x;d) 


Kiang PDF , formula ([8]) when d = 


2 , 


2142 


H{x;d) 


Kiang PDF , formula ([8]) when d ~ 


0.61 , 


43 


h{v) 


product two gamma, formula (jlSp 




40 



4. Non-Poissonian seeds 

We now explore the case in which the seeds of the Voronoi Diagrams are 
distributed in a correlated way with respect to the center of the box. The 
correlated seeds are generated in polar coordinates , p and 9. The radius p , 
the distance from the center of the box , is generated according to the PDF 
which is the product of two GV2 , see formula once the scale parameter b 
is introduced , 

hiv,b) = . (23) 

Random numbers of the polar angle 9 , in degrees, can be generated from those 
of the unit rectangular variate R using the relationship 

61 - 360 X i? . (24) 

A typical example of Voronoi Diagrams generated by such seeds is reported 
in Figure [5] and Table ([5]) reports the oi three different fits. From a careful 
analysis of Table ^ it is possible to conclude that the product of two random 
variables, produces a better fit with respect to the Kiang function with fixed or 
variable d. 

Another example is represented by the ratio of two GV2 , formula ((2T|) , 
that when the scale parameter b is introduced has PDF 

6 vb"^ 

r{v,b) = 5 . (25) 

b-^{b + vf [h^ -^2vb + v-^) 

An example of Voronoi Diagrams generated by seeds that follow the quotient 
of two GV2 is reported in Figure [7] and Table [6] reports the of three different 
fits. 

From a careful analysis of Table [H] it is possible to conclude that the quotient 
of two GV2 , also in this case produces a better fit with respect to the Kiang 
function with fixed or variable d. 

A third example is represented by the Normal ( Gaussian ) distribution 
which has PDF 

1 

^^""'"^^ " a(2^)i/2 '^'^P(~2^) -oo<x<oo . (26) 
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Figure 6: The Voronoi Diagram in 2D with Non-Poissonian seeds generated according to 
the product of two GV2 , see formula Ij23|l . The selected region comprises 42 random seeds 
marked by a point and 6 = 25. 



Table 6: Normalized area distribution in 2D of 117 correlated seeds generated according to 
formula I I25I I. The of data fit when the number of classes is 40 for three PDF 



PDF 


name 






H{x; d) 


Kiang PDF , formula ^ 


when d ~ 2 , 


5560 


H{x; d) 


Kiang PDF , formula dH) 


when d = I , 


87 


r{v) 


ratio two gamma, formula 


m 


79 
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Figure 7: The Voronoi-diagram in 2D when the seeds are generated according to the quotient 
of two GV2 , see formula l|25|l . The selected region comprises 117 seeds and b = 25. 
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Tabic 7: Normalized area distribution in 2D of 145 correlated seeds generated according to 
formula l|27| l. side=l, cr = 0.5. The of data fit when the number of classes is 40 for three 
PDF 

PDF name 

H{x] d) Kiang PDF , formula © when d^ 2 , 2070 
H{x; d) Kiang PDF , formula dH]) when d = 0.88 , 84 
HN{x;a) Half Normal{Gaussian) PDF, formula ^ 168 
hn{v] b) Product two Gaussian PDF, formula (|29|) 78 



When only positive values of x are considered PDF ([26]) transforms in 

^^^'^=^^ = ^(2^"''^^"^) 0<a:<oo . (27) 

We now consider the product of two normal random variables X — N{0, 1) 
and Y = N(0, 1). The PDF of V = XY is ^ 

M, f Ko(v * siqnum(v))/TT —oo<v<0 
^ S ) ) \M c (28) 

I Kq[v * signum(v ))/Tr u < v < oo . 

This PDF has a pole at u = and when the scale parameter b is introduced and 
only positive values are considered it transforms in 

hn{v; b) = 2Ko{^)/t: Q < v < oo . (29) 

An example of Voronoi Diagrams generated by seeds that follow the Half 
Normal (Gaussian) distribution (formula ([27ll ) is reported in Figure [8] and Ta- 
ble [7] reports the of four different fits. 

From a careful analysis of Table [7] it is possible to conclude that the product 
of two Gaussian PDFs produces a better fit with respect to the Half Normal 
(Gaussian) distribution and the Kiang function. 



5. Aggregation of molecules 

The analysis of the molecular dynamics through the Voronoi polyhedron in 
3D or the Voronoi polygons in 2D is becoming a standard procedure. What is 
very interesting is the surface aggregation of a water-methanol mixture at 298 K, 
see [131 • The aggregation of the methanol and water molecules at the surface 
of their mixture is analyzed in the light of the Voronoi diagrams conside ring 
water and methanol molecules together , see Figure 11 in [30]. Figure 12 in [30| 
reports an instantaneous snapshot of the surface layer of a system containing 5 % 
methanol molecules (top view), as taken out from an equilibrium configuration. 
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Figure 9: Superposed on the V — P network of Figure [2] are present 131 molecules of water 
(red small balls) and 54 molecules of methanol (green great balls). 
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Figure 10: Superposed on the network of Figure |6] generated by non Poissonian seeds are 
present 131 molecules of water (red small balls) and 54 molecules of methanol (green great 
balls). 



In order to simulate such an aggregate Figure [2] reports the 2D V — P and 
Figure [9] the theoretical displacement of the molecules on such a network. 

The same aggregate is simulated by non Poissonian seeds and Figure [6] re- 
ports the Voronoi Diagram while Figure \W\ the theoretical displacement of the 
molecules on such correlated network. 



6. Conclusions 

The size distribution of — P in ID has an exact analytical result that is 
a GV2. Starting from this analytical result we derived the sum , the product 
and the quotient of two GV2. The area distribution of F — P in 2D is well 
represented by the sum rather than the product of and quotient of two GV2 , 
see Table [T] A careful comparison between the Kiang function, formula ([5]) , 
and a new function suggested by Ferenc-Neda , see formula (|13p , has been done 
, see Table [T] and Table [S In presence of non-Poissonian seeds that present 
a symmetry around the center the situation is inverted and the product of 
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two GV2 , formula ^5]) . or the ratio of two GV2 , formula (PT|) , produce a 
better fit of the area distribution of Voronoi polygons in 2D with respect to 
the Kiang function, formula ([5]) , see Table [5] and Table [HI A third test made 
on Normal (Gaussian) correlated seeds is in agreement with the conjecture that 
the area of the 2D Voronoi diagrams follows the distribution of the seeds. The 
developed algorithm allows us to simulate some well studied aggregation such 
as the molecules of the water-methanol mixture. The question of whether if 
the area of the voids between molecules follows the V — P area distribution as 
represented by the Kiang function , see equation ([5]) or the gamma of Ferenc- 
Neda , see equation (fT5|) is open to future efforts. 
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